ON COMPUTING HIGHER-ORDER ALEXANDER MODULES OF KNOTS 



PETER D. HORN 

Abstract. Cochran defined the nth-order integral Alexander module of a knot in the three sphere as 
the first homology group of the knot's (n-fl)th-iterated abelian cover. The case n = gives the classical 
Alexander module (and polynomial). After a localization, one can get a finitely presented module over 
a principal ideal domain, from which one can extract a higher-order Alexander polynomial. We present 
an algorithm to compute the first-order Alexander module for any knot. Included in this algorithm is 
a solution to the word problem in finitely presented Z[Z] -modules. 



1. Introduction 

Given a diagram for a knot K in S 3 , there is a well-known algorithm for producing the Wirtinger 
presentation of tti(S 3 \ K) [Rol76, Section 3.D]. While the knot group is a complete knot invariant 
(see [Wal78]), using these presentationsis an impractical tool for distinguishing knots. Since any knot's 
exterior has Hi(S 3 \ K) = Z, abelianizing the fundamental group is not at all a useful invariant for 
distinguising knots. A logical next step would be to find an intermediate quotient, Q(K), that is more 
discerning than Hi and more tractable than tt\: 

7n(5 3 \ K) -» Q(K) -» Ht(S 3 \ K) 

Or better yet, one may hope to find a module M(K) on which the group Q(K)-acts. With this strategy 
in mind, Cochran defined the higher-order Alexander modules of a knot [Coc04]. The derived series 
of a group G is defined by G^ — G and G'" +1 ' = [G^ n \ G^] , where the square brackets denote 
commutators. Throughout this paper we are interested in knot groups G — tti(S 3 \ K). Put simply in 
algebraic terms, the nth-order integral Alexander module of K is G^ n+1 ^ /C?( n + 2 ) as a rig htZ [G/G {n+1) ]- 
module (the action is conjugation in the group). For n — 0, this is the classical Alexander module, 
which has a topological interpretation as the first homology of the universal abelian cover of the knot 
complement. The higher-order Alexander modules we consider here are different from those defined 
before in Cochran-Orr-Teichner [COT03, Sections 2 and 3]; the Cochran-Orr-Teichner modules are 
defined using a localization stemming from a coefficient system on a 4-manifold in which the knot K is 
slice. The modules considered here only depend on the knot K . 

Let A\(K) denote Cochran's nth-order integral Alexander module. In [Coc04], Cochran constructs 
a localization of the coefficient ring and considers the nth-order (localized) Alexander module A n (K). 
A precise definition will follow in Section 2. As in the classical case, one can extract an integer-valued 
invariant 5 n (K), which is the degree of an "nth-order Alexander polynomial." These degrees give lower 
bounds for the genus of a knot, which is the minimal genus of all orientable surfaces in S 3 whose boundary 
is the given knot. 

Theorem 1 (Theorems 7.1 and Corollary 9.2 of [Coc04]). If K is a knot in S 3 with non-trivial Alexander 
polynomial, then Sq(K) < 5i(K) + 1 < • • • 5 n (K) + 1 < • • • < 2genus(if). Furthermore, given any n > 0, 
there exists a knot where 5 n (K) < 5 n +i(K). 

Thus, the higher S n give better lower bounds for the genus of a knot. For the '<' part of the theorem, 
Cochran constructs knots by a satellite operation called 'infection.' Many positive results in the so- 
called 'higher-order knot theory' have been proven using satellite operations, but satellite knots are not 
generic. The higher-order degrees 5 n are the most alluring of all the higher-order knot invariants for 
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two reasons: they are denned for all knots, and they 'should be' algorithmically computable. Since 
the S n were defined, several people have made successful computations by hand [LM08], [HollO]. In a 
related line of inquiry, Harvey successfully computed degrees of higher-order Alexander polynomials of 
3-manifolds [Har05, Sections 6 and 8]. Such calculations are tedious but always seem to work out with 
enough perseverance. 

The purpose of this paper is to describe a practical algorithm to compute Si(K) for any knot K, 
answering question 13 of [FV11]. Using this algorithm, we compute Si for many low crossing knots. 
All knots with ten crossings or fewer have So — 2 genus (if), and so we know Si by Theorem 1. It is 
easily seen from the definitions in Section 2 that if the classical Alexander polynomial of K is trivial, 
then all the S n {K) = 0. The remaining knots with eleven crossings for which < So < 2 genus are 
11„4 5 , ll n 67) Hn73i Hn97i an d H n i52- We exhibit the utility of the algorithm to compute <5i(ll„67) and 

<5l(12 n 293)- 

Theorem 2. <5i(ll n 67) = 3 and <5i(12„293) = 3. 

The classical Alexander polynomial of 11„67 is 2 — 5t + 2t 2 , and its genus is 2. This is the simplest 
known example of a knot where Si is a better bound on the genus than So. (The only simpler possible 
example is 11„45). Since 3 + 1 = 2-2 (as in Theorem 1), we conclude that S n (ll n Q^) = 3 for all n > 1. 
We remark that the classical Alexander module of 11„67 is cyclic; it matches that of the stevedore knot, 
6i, which has genus one. The calculation of <5i(lln67) took 11324 seconds on a machine with an Intel i5 
2.5 GHz quad-core processor. 

The knot 12„293 also has genus 2, and its classical Alexander polynomial is 2 — 3t + 2t 2 . 

For each of 11„67 and 12 n2 93, there exists another knot in the tables with isomorphic knot Floer 
homology; the computation was verified by Gridlink [Cul]. These other knots are ll n 97 and 12 n 5ig, 
respectively. We were unable to find Wirtinger presentations for these knots that would allow our 
program to compute Si in a reasonable amount of time, which leaves the question: does knot Floer 
homology detect Si? 

In their survey on twisted Alexander polynomials [FV11], Friedl and Vidussi ask whether the higher- 
order Alexander polynomials detect mutation, and as of the writing of this paper, their question remains 
open. 

2. Definitions 

The definitions of higher-order Alexander modules are due to Cochran [Coc04] . 

Let if be a knot in S 3 and G = tti(S 3 \ K). Let Y n = G/G<- n+1 \ where G( n+1 ' is the (n + l)th term 
of the derived series of G. There is a rig ht action of T„ on G {n+1 ^ /G ( ™+ 2) by conjugation: 

g^h^h- 1 gh, for all h € r„, g G G (n+1) /G {n+2) 

Let M n denote the (n + l)th iterated (maximal) abelian cover of S 3 \ K, so that tti(M„) ^ G(™ +1 ). 
Then Hi(M n ;Z) = G^ n+1 ^ / G^ n+2 \ and this homology group has a right action by r„ (viewed as either 
deck translations or conjugation). 

Definition 1. The nth-order integral Alexander module of K, denoted A^(K), is the right Zr n -module 

Al{K) := g(" +1 Vg ( " +2) = Hi(S 3 \ K;ZT n ) S Hi(M n ;Z) 

The group r„ is poly- (torsion-free) abelian (PTFA), i.e. it admits a scries of normal subgroups 
{e} < Ak < Ak_i <•••<] A = r„ so that each Ai/A i+ i is torsion-free abelian. To see this, take 
Ai = G^/G {n+1 ^ and recall that Ai/A l+ i = G^/G < - 1+1 ^ is torsion-free abelian [Str74]. We will use 
the fact that Zr„ imbeds in its classical quotient field (in general 2T„ is a noncommutative ring, so its 
quotient field is a skew field): 

Proposition 1 (Proposition 3.2 of [Coc04], see also the Proposition of [Lew72]). IfT is PTFA, then QT 
(and hence XT) is a right Ore domain, i.e. Qr imbeds in its classical right ring of quotients JC, which is 
a skew field. 
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Let K. n denote the quotient field of ZT n . There is an important intermediate ring R n (Qr„ C R n C 
/C„) which will be convenient for our use. For n > 1, the kernel of tt : G/G^ -> G/G^ = Z is just 
After choosing a splitting 1 i — y [i of 7r, we get an isomorphism 

r n = G/G (,i+1) S G (1) /G (n+1) xi Z 

Any element of r„ has a unique expression as fi l g, where 5 € G^/G^ 1 ). There is also a unique 
expression by writing powers of fj, on the right: /j, 1 g = g fi l , where the ~ denotes the action on g by 
conjugating i times by /1. Thus QT n is canonically isomorphic (after choosing the splitting) to the skew 
Laurent polynomial ring (Q [G^ 1 ' / G^ n+1 ']) [* ±:l ]< In this skew polynomial ring, the coefficients are not 
commutative (unless n = 1), and the coefficients do not commute with the variable (due to the semidirect 
product structure). 

Let K„ denote the classical right (skew) field of quotients for Q [G^ /G^" +1 - ) ] , and let 

Rn = Zr„ (Z [G« /G(" +1 '] \ {0})~ . After choosing a splitting f 1 y /i, there is a canonical isomorphism 

Rn K„ [t^ 1 ]. Since R n is a localization of Zr„ wherein some but not all elements are inverted, we 
have that Zr„ C Rn C K, n . 

Definition 2. The nth-order (localized) Alexander module of K is A n (K) — Hi(S 3 \ K;R n ). 

We now summarize several useful facts about R n and IC n : 

Theorem 3 (Section 4 of [Coc04]). Suppose we have chosen a meridian (i.e. splitting) 1 4/1 so that 
R n S K n [t ±l ]. Then 

(1) R n is a right (and left) PID, 

(2) K„ [i ±:L ] has a well-defined degree function and a Euclidean algorithm, 

(3) A n {K) is a finitely-generated torsion module over R n , 

(4) R n is a flat left ZT n -module, so that A n (K) = A„(K) ®zr n Rn, 

(5) K n is a flat R n module so A n {K) ®R n K n — Hi(S 3 \ K; JC n ) 

The structure theorem for finitely-generated modules over commutative PIDs generalizes to the non- 
commutative case [Coh85]: if M is a finitely generated torsion right R module, where R is a PID, then 
M = R/ei R © • • • © R/ek R where is a total divisor of e^+i (moreover, this determines the up to 
similarity) . 

If one views the ei as elements of K„ [t ±:L ] , one could define the higher-order Alexander polynomials 
of K as the product of the e^. However, this is not a well-defined invariant of K, as the isomorphism 
R n = K„ [t ] depends on a choice of meridian. However, the degree of this polynomial is a knot 
invariant. 

Definition 3 (Definition 5.3 of [Coc04]). For any knot K and any n > 0, the degree of the nth order 
Alexander polynomial, denoted 5 n (K), has several equivalent definitions: 

(1) the sum of the degrees of the ei € R n = K„ [t ] , 

(2) the rank of A n (K) as a module over K„, 

(3) the rank of G {n+1 ^ /G {n+2 ^ <g) zr „ Rn as a module over Z [GW/G ( " +1) ] , or 

(4) the rank of G {n+1 *> /G {n+2 *> as a module over Z [G^ /G < - n+1 ^] 

Example 1. If n = 0, r = G/G^ Z. In this case K = Q, and R ^ Q^ 1 ]. Then A%(K) 
is the classical Alexander module G^/G^ over Z[t ±l ], and Aq(K) is the rational Alexander module 
Ao(K) = GW/GW ® Z [f±i] Q[* ±x ]- 

Example 2. Consider the case n = 1, which is the main focus of this paper. The group ri = g/g< 2 > 
fits into an exact sequence 

o^g (1) /g (2) ->r a ^z^o 

After choosing a preferred meridian, which will be identified with t, we see that any element of Ti can 
be written uniquely as t % g, where g £ G^/G^. In this case, Ki is the quotient field of Z [GW/G (2) ] , 
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which is a commutative field! Addition and multiplication in this field work - at least on a symbolic 
level - as they do in Q, but one must be especially careful not to divide by zero. Section 4 addresses this 
concern. While Ki is commutative, the coefficients do not commute with the variable in i?i = Ki[£ ±:l ]. 
For example, tg = \ig\i~ l t in the polynomial ring. 

2.1. Outline of the algorithm to compute 5\. To compute 5i(K), we will take a presentation matrix 
M for A\(K) as a module over K^^ 1 ] and make it upper-triangular. The sum of the degrees of the 
diagonal entries is 8\ (K) . Note that we do not need the divisibility condition on the diagonal entries as 
in the structure theorem (cf. Lemma 1). The triangularization process requires only row operations (and 
some column swaps), which can be achieved by left-multiplication by elementary matrices over Ki[i ±:L ]. 
The algorithm to upper-triangularize a matrix over ] is the same as over Q[i ]. First, multiply 

each row of M by a suitable power of t to make all entries honest polynomials. Pick the smallest degree 
entry and move it to the top-left entry by row and column swaps. Use row operations to clear the first 
column; it may be necessary to perform a row swap if the degree of an entry in the k, 1-position (where 
k > 1) falls below the degree of the 1, 1-entry. By doing these row operations, the matrix will eventually 
have a non-zero entry in 1, 1 and all zeros below that; this follows from items 2 and 3 of Theorem 3. 
The algorithm continues by moving on to the submatrix gotten by removing the first row and column. 

The difficulty with implementing this algorithm lies in the coefficient field Kj. , which consists of formal 
quotients of formal Z-linear combinations of elements in the classical Alexander module, G^'/G^ 2 ', the 
group operation for which we will denote additively in this paragraph. While this is a commutative field, 
arithmetic is lengthy. For example, the sum 

ma + nb pd mq(a + e) + mr(a + /) + nq(b + e) + nr(n + /) + op(c + d) 

oc qe + rf oq(c + e) + or(e + /) 

seems to have 5 terms in its numerator and two terms in its denominator. But do any of the terms 
cancel? Is the numerator zero? Logically, the 'collapsing problem' in the group ring ZH is equivalent 
to the word problem in the group H. However in practice, one would prefer to have a 'normal form' for 
group elements in H to speed up the ring operations in ZH. This will be discussed further in Section 4. 

While Ki is a commutative field, its arithmetic seems to push the boundaries of the abilities of modern 
computers. It seems to the author that an implementation to compute 5 n for n > 2 is hopeless with the 
current technology, since the fields K n with n > 2 are noncommutative. 

3. Previous ad hoc computations 

For knots with very few crossings (five or fewer), one can calculate 5% by hand. Such a calculation 
is straightforward for two reasons: the size of the matrix you must work with is manageable, and the 
coefficient field Kj is simple, due to the simple structure of Aq for low crossing knots. We present two 
examples to illustrate both the methods involved in the general algorithm and why previous 5\ compu- 
tations have failed to give a general algorithm. Compare the strategy in the next example with [HollO, 
Example 4.2]. 

Example 3. Let T denote the trefoil knot. Using the standard torus knot diagram, one can calculate 
the Wirtinger presentation 

G := 7Ti (S 3 \ T) = (xi,x 2 ,x 3 \ x i x 3 x^ 1 x 2 1 , xix 3 x^ 1 x^ 1 ) 

For reasons that will soon become clear, change the generating set to a := X\, y\ := x 2 x^ 1 , y 2 := x^x^ 1 . 
Note that under the abelianization map, a maps to a generator of Z while y\ and yi map to the identity. 
We now have the presentation 

vri (S 3 \T) = (a,y 1 ,y 2 \ ay 2 a~ 1 yY 1 , ay 2 yi 1 a^ 1 y 2 ' 1 ) 

We construct the matrix of Fox derivatives [Fox53] whose i, j'-entry is the derivative of the ith relation 
with respect to the jth generator; the Fox calculus is defined by 

• jB = <%> for each hi 
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• S^- = 0, for each i 

Ogi ' 

• = + w J^~7 f° r each i and for each u,v in the group 

where the are the generators of the group. It follows that = -u- 1 -^ for all u in the group. The 
Fox matrix for the second presentation is 

/ 1 - ay 2 a~ 1 -1 a \ 

The entries of this matrix lie in the group ring Z7ri . Upon writing the group elements in any quotient 
group G, one obtains a presentation matrix for Hi (S 3 \ T, xq\ ZG), the first homology relative to a 
basepoint of the covering space whose group of covering transformations is G. 
Let us write the Fox matrix 'over Z' by abelianizing the group elements: 

(0 -1 t\ 

The first column has all zeros by the choice of generators a, 2/1,2/2- Recall that the Fox matrix is a 
presentation matrix for the first homology of a covering space relative to a basepoint; this, in effect, 
adds a free generator to Hi (S 13 \ T;Z \t ]), since in the long exact sequence for the pair (S 3 \ T,xq), 
we have 



i/ 1 (^ 3 \T;Z[t ±1 ]) i/ 1 (5 3 \T,x ;Z[i ±1 ]) H (x ; Z [i* 1 ] ) 



torsion free © torsion Z [i* 1 ] 

By exactness, the free part of Hi (S 3 \T,xo;Z [t^] ) has rank one. 

The free generator in this case is the one corresponding to the first column, and upon eliminating 
that column, we end up with a presentation for Hi (S 3 \ T; Z It 1 ) , i.e. the classical Alexander module 
of T: 




By the definition of 2/1,2/2, V\-,Vi may be viewed as generators of the classical Alexander module (they 
lie in the commutator subgroup), and the first and second columns of the above matrix correspond to 
these generators. In other words, Aq(T) is generated as a Z [t^ 1 ] -module by y\ and 2/2 with the relations 
— 2/1 + t * 2/2 = and — t * 2/1 + (t — 1) * 2/2 = 0. One can actually 'upper-triangularize' this presentation 
matrix (using only row operations, which amounts to adding multiples of the relations together) to: 

(0 t 2 -t + l) 

Thus 2/1 = t * 2/2 and (i 2 — t + 1) * 2/2 =0 in the Alexander module Aq(T). This description of Aq(T) 
will allow us to understand the coefficients in Ki. 

To compute Si(T), we need a presentation matrix for Ai(T) = Hi (5 3 \ T;7Zi). Since ZTi C TZi, 
and Ti is a quotient of 7Ti, such a presentation matrix can be derived from the Fox matrix above. First, 
reduce the elements in the Fox matrix according to the quotient map 7Ti (£ 3 \ T) -» G/G^ = T±. Recall 
from Section 2 that any element in Ti has a unique expression (using the semidircct product) once we 
have chosen a meridian. We have already chosen a distinguished meridian: a. Let us write the Fox 
matrix with entries written in G/G^ G«/G (2) * Z, with t generating Z (we use a when it conjugates 



6 



PETER D. HORN 



an element of G^ /G^ 2 \ and t when there is no conjugation): 

1 — ay 2 a~ 1 —1 t 

1 - ay 2 a^ 1 ay^ 1 a~ 1 -ty 2 y^ x t—l 

Again, this matrix presents Hi (S* 3 \ T, xo; ZFi). Since Ri is a flat ZTi module, so this matrix also 
presents Hi (S* 3 \ T, xq; Ri) since ZFi To get a presentation matrix for Hi (S 3 \ T; i?i), we 

must eliminate one column as we did above for Z [t ±x ] . Let denote the ith relation and Xj the jth 
generator in the presentation for G. By the fundamental formula for Fox derivatives [Fox53, eq 2.3], 

This implies that the jth column of the Fox matrix is a i?i-linear combination of the other columns 
as long as x,- € G (1 ) and Xj ^ 1 e ITi. As long as we can find a generator y, k 6 G^ 1 - 1 which is not 
in we can (perform column operations according to Equation 1 until the column corresponding 

to yi has all zeros and) delete the column corresponding to yi. This results in a presentation matrix 
for Hi(S 3 \ T;Ri), since rank^ffi (S 3 \ T, x ; Ri) = 1, which we will justify now. The existence 
of the column of zeroes establishes that rank^iJi (S* 3 \ T, xo',Ri) > 0. One may choose a splitting 
Hi (S 3 \ T, xo', Ri) — Rf © torsion and consider the map (Ri) d — > Hq(xo', Ri) — Ri induced by the right 
most map in the following portion of the long exact sequence for the pair (S 3 \ T, Xq): 

H 1 (S 3 \T;R 1 ) Hi(S 3 \T,x ;Ri) H (x ;Ri) 



torsion free © torsion i?i 



The map (Ri) d — > Ri is injective by the exactness of the sequence above. That d = 1 now follows from 
the rank- nullity theorem over a PID. 

The upper-left triangular presentation matrix for Aq(T) guarantees that y 2 ^ 1 G G' 1 -* /G^K Thus, 
we remove the last column of the Fox matrix over ZTi (which corresponds to y 2 ) to obtain a square 
matrix, which we will call the metabelian Fox matrix: 



(2) 



1 - ay 2 a 1 -1 
1 - ay 2 a~ 1 ayi 1 a~ 1 -ty 2 y^ 1 / 

The reader may have noticed that the second column may also be deleted (the element yi is nontrivial 
in GW /G' 2 )). We chose the last column because fewer row operations are required to upper-triangularize 
matrix 2. 

Since ZTi c TZi = Ki [t^ 1 ], matrix 2 is a presentation matrix for Ai(T), as discussed above. We 
choose in our algorithm to use only row operations at this point because they are achieved by multi- 
plications on the left by elementary matrices over Ki Column operations are equally valid but 
require the extra mental effort of remembering to do multiplication on the right in the noncommutative 
ring Ki [t ±x ]! 

We add —t y 2 yi X times the first row to the second and switch columns: 

-1 1 - ay 2 a~ 1 

t (ay 2 a~ 1 y 2 y^ 1 - y 2 yf 1 ) + 1 ~ ay 2 a~ 1 ay^ l a~ 

At this point, adding the degrees of the diagonal polynomials gives #i(T) by Lemma 1, but we must 
be careful. The —1 has degree zero. Linear- looking polynomials such as the bottom-right entry do not 
necessarily have degree one; for example, tx + y will have degree 1 if and only if y ^ and i^0£ Ki. 
Recall that Ki is the quotient field of Z [G« /G^] . We can determine whether the coefficeiits are zero 
using the presentation matrix for Aq(T). The i-coefficient ay 2 a~ 1 y 2 y^ 1 — y 2 y-[ 1 is zero if and only if 
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ay 2 a~ 1 ayi 1 a~ 1 = 1 in G^/G^ 2 \ This equality in the group can be translated into an equality in the 
Z [t^] -module Aq(T): t*y 2 — t*yi =0? Since y\ ^ y 2 in Aq(T), the ^-coefficient in the linear-looking 
polynomial is nonzero. We leave it to the reader to verify that the constant term is also nonzero. 
We compute 5 X (T) = deg(-l) + deg(t {ay 2 a~ 1 y 2 y^ 1 - y^ 1 ) + 1 ~ ay 2 a~ l ay^[ l a~ l ) = + 1 = 1. 

Few choices were made in the computation of S\(T) (choice of preferred meridian, which column 
to delete). These choices did not affect the end result, and one can easily program a computer how 
to make these choices. The burden of the computation lies in putting the metabelian Fox matrix in 
upper-triangular form. In general, many row (and perhaps column operations) are needed, and the 
Ki-coefficients quickly become unruly. While Example 3 was fairly simple, an outline for the general 
algorithm to compute 5i can be gleaned: 

Algorithm 1. 

(1) Input: a diagram for a knot K. 

(2) Compute a Wirtinger presentation for tti(S 3 \ K) and change the generators. 

(3) Compute the general Fox matrix and abelian Fox matrix. 

(4) Be able to recognize in the ring ZFi so row operations over Ki \t ] are possible. 

(5) Decide which column(s) of the Fox matrix can be deleted, and delete one of them to arrive at 
the metabelian Fox Matrix. 

(6) Put the metabelian Fox matrix in upper-triangular form using valid row operations over Ki [t ±:L ] ■ 

(7) Compute and add the degrees of the polynomials on the diagonal, giving 8\{K). 

Proof. We will describe an algorithm for each step and a lemma pertaining to steps 6 and 7. Recall 
that Definition 3 item 1 requires a diagonal matrix that presents A\ (K) wherein the diagonal entries are 
total divisors of the subsequent diagonals. The divisibility criterion is immaterial in the computation 
of the degrees of the diagonal polynomials. Lemma 1 guarantees that an upper-triangular, rather than 
diagonal, form of the presentation matrix suffices, which will save considerable time by avoiding many 
tedious row and column operations. 

Steps 1, 2 and 3 are straightfoward. We use in our implementation a script written by an REU group 
at Columbia University that converts a knot diagram (PD Code) to a Wirtinger presentation [REU11]. 
Step 4 is the most difficult. One needs to understand Z [G^/G^] well enough to know when an 
element is zero or nonzero; outside of that, addition, multiplication and inversion in Ki behave like the 
same operations in Q. The proofs of Theorem 4 given in Sections 4.1, 4.2, and 4.3 settle step 4. Step 
5 is then simple as in Example 3: the column corresponding to one of the generators yi € G^ 1 ' may 
be deleted if and only if yi ^ in /G^ 2 \ Step 6 is settled by an adaptation to Ki [t^ 1 ] of the 
algorithm to compute the canonical form of a module over a PID (cf. Section 2.1), although for the 
purpose of computing <5i , one only needs upper triangular and does not need the divisibility criterion as 
discussed above. Finally, step 7 is straightforward (we use the degree function for Laurent polynomials, 
so deg(t + t- 2 ) = 3 and deg(< 2 + t) = 1). 

□ 

Several groups of undergraduates have attempted to implement an algorithm to compute <5i , and the 
sticking point tends to be step 4. A strategy that works for many low crossing knots is to compute 
the abelian Fox matrix, which is a presentation for Aq(K) = G^ /G^ 2 \ and put it in upper triangular 
form. This works for the trefoil in Example 3. The problem is that 'upper-triangularizing the abelian 
Fox matrix' is impossible in general, as we see now. 

Example 4. Let K denote the knot in Figure 1. One may compute via Fox calculus or from a Seifert 
surface a presentation for Aq(K) to be 

_ /3 - 3t 2-t 
\l-2t 3 - 3i 

The Alexander polynomial of K is It 2 — 13t + 7, an irreducible polynomial over Z [t^ 1 ] ■ 
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P can be put in upper-triangular form over Q [i 1 , which is a PID, but this is impossible over Z [t ] , 
for suppose 




where ~ denotes a sequence of invertible row and column operations over Z[i ±r |. We may assume 
without loss of generality that a = 1, since 7i 2 — 13£ + 7 is irreducible. 

Recall that the first elementary ideal of a module with an r x r presentation matrix is the ideal £ i 
generated by all r — 1 x r — 1 minors of the matrix. For a 2 x 2 matrix such as P, E\ is the ideal generated by 
all of its entries. This ideal is invariant under the operations ~. For our P, E\ — (2 — i, 1 — 2t), since 3 — 3t 
is the sum of the two off-diagonal entries. If P were upper-triangularizable as above, then £ i = Z [£ ±:L ] ■ 
Then we could write some power of t as t m — p(t) (t — 2) + q(t) (2t — 1) for some p(t), q(t) £ Z[t]. Plugging 
in t = 2, we see that 2 m = q{2) * 3, but 2 m is not divisible by 3, a contradiction. 




Figure 1. A knot with non- upper-triangularizable abelian Fox matrix 



We now present the lemma which is crucial in the proof of Algorithm 1. 

Lemma 1. Suppose P is a square presentation matrix for A n (K), or more generally for any torsion 
K n [t* 1 ] -module. If P is upper-triangular with diagonal entries c?i, . . . ,dk, then 

k 

^dcgdi = rank Kn .4„(.?0 = S n (K) 
t=i 

Proof. The claim follows essentially from the Euclidean algorithm. Row (respectively, column) opera- 
tions are realized by left (respectively, right) multiplication by elementary matrices. Thus, to reduce an 
entry c in the matrix by another entry b of smaller degree using a row operation, we need to use the 
'left quotient' c = qb + r. Right quotients will be used for column operations. Since K„ [t^] is a PID, 
one can transform P to a diagonal matrix using these operations. Furthermore, it suffices to prove the 
lemma for 2 x 2 matrices. For, suppose the lemma has been established for 2 x 2 matrices. Given 

(a b c 
Ode 

/ 

we may change P to 

fa' c 

P ~* d' e' 

\0 / 
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by only disturbing the first two rows and columns. This may in turn be changed to 




by extending the operations on the submatrix 

a 1 c 
/. 

to the entire 3x3 matrix. One continues until P has been completely diagonalized. One continues by 
induction to prove the lemma for matrices of arbitrary size. 
It remains to establish the 2x2 case. Suppose 



P = 



a b 
c 



Given an element x e K„ [i* 1 ], we denote its degree by using uppercase symbols, i.e. X = dcg.x. 
Returning our focus to P, we may perfome some row and column operations to guarantee that B < A 
and B < C. The ring K n [i ±:L ] has the Euclidean algorithm. We compute the quotients qi, ■ ■ ■ q n +2 and 
remainders n . . . r n+ \: 

c = qib + rx 
b = q2T\ + r 2 
r\ = q 3 r 2 + r 3 



r n -i = q n +ir n + r n +i 
r n = q n +2r n +i + 

Let Oi for i = 1, . . . , n + 2 denote the row operation which adds — qt times row | + §(— 1)* to the 
other row. The first operation is 







( a 








\-qia 





Note that 0\ takes P out of the set of upper-triangular matrices, but the full sequence of the O, almosts 
takes P back to upper-triangular. To see this, define a recursion relation by /_i = 0, fa = a, and 
fk = fk-2 — qkfk-i- We see immediately that the degrees satisfy Fi = Qi + for i > 1, since Qi > 0, 
and we see that 

a b\ o, ffo b\ Ch f.f2 r 2 \ (fc r 2 \ cu ( ,f 4 r 4 
c) \h rj {.h rj {.h r 3 J {f 3 r 3 

Since r n+2 = 0, this sequence of operations ends at either 

fn+i r n+1 \ ot //„ +2 

fn+2 J \fn+l r n+ i 

depending on the parity of n. Let S denote the operation which either swaps the two columns or swaps 
the columns and swaps the rows, so that 

fn+l fn+1 
fn+2 



S(O n+2 ■ ■ ■ OxP) 
is upper-triangular. Perform a column operation to arrive at 



pi _ I r n +i fn+l 
' f n+2 



10 



PETER D. HORN 



where F' n+1 < R n +i < B. Note that R n+ \ + F n+2 = R n +i + Qn+2 + F n +i = R n + Fn+i by the definitions 
of the Ti, qi, and /j. We apply this observation repeatedly to conclude that R n +i + F n+ 2 = Ri + F2, but 
this is in turn i?i + Q2 + Q\A = B + Qi + A = C + A, which is the sum of the diagonal degrees of P. 

Thus, we have described a procedure that takes P = [ n to P' = " , such that A + C = 

A 1 + C and B > B' . After repeatedly aplying this procedure, we end at a diagonal matrix (b 1 = 0) 
which presents the same module as P and has the same diagonal degree sum. □ 

4. The word problem in finitely presented Z[Z]-modules 

Recall that Ki is the (commutative) quotient field of the ring Z [G^'/G' 2 '] . In doing the Ki [t ±1 ]-row 
operations to put the metabelian Fox matrix in upper-triangular form, it will be necessary to compute 
(left) polynomial quotients in the ring Ki [t^ 1 ] ■ Thus, one needs to invert elements of Ki and in 
particular, needs to know when an element in Z [G^-'/G'- 2 ''] is zero so as not to divide by zero in Ki! 
This is the difficulty of step 4 of Algorithm 1. 

Let H be any group, and consider its group ring ZH. Given an element p = n\h\ + • • • + n^hk of 
Zif , one may determine whether or not p — by the algorithm: for 2 < i < k, if hi = hj for some j < i, 
combine the ith and jth terms by adding rij + Uj. After running this 0{k 2 ) algorithm, p — if and only 
if all remaining m — 0. 

Suppose we have a presentation of the group H with generators Xi and relators rj (neither indexing 
set must be finite), and let h be an arbitrary element of H, written as a word in the generators. The 
problem of deciding whether h is equal to the identity element (represented by the empty word) is known 
as the word problem for this presentation of H . We say this presentation of H has solvable word problem 
if there is an algorithm taking any element h = h(xi) of H and deciding whether ft, = 1 or not. We 
refer the reader to Stillwell's enlightening survey on the word problem [Sti82] . While much effort has 
been made to understand the word problem for finitely presented groups, the case of infinitely presented 
groups has received less attention. 

We are interested in the group H = G (1) /G (2) , where G = n^S 3 \K). In this case H is an abelian 
group, and H is finitely generated if and only if the (classical) Alexander polynomial Art (t) is monic. 
The word problem for finitely generated abelian groups is solvable by the structure theorem for finitely 
generated Z-modules, but sadly most knots do not have monic Alexander polynomial. Recall that H 
has the structure of a Z[Z]-module. Even though Z[Z] is not a PID (and so its modules have no a priori 
structure theorem), the word problem for such modules is solvable. 

Theorem 4. Let H be an abelian group. Suppose that H has finite presentation P as a module over 
Z[Z]. Then there exists a group presentation P' for H that has solvable word problem. 

Note that P is not a presentation of H as a group, though a group presentation can be obtained 
from P. For each Z[Z]-module generator x, one considers the Z-many generators t l * x. Similarly, each 
Z[Z]-module relation gives rise to Z-many group relations. These generators and relations give a group 
presentation P' for H . 

The author could not find this result in the literature and expects it to be 'folklore.' The special case 
when H is the classical Alexander module of a knot in S 3 is well-known to knot theorists; see Sections 4.1 
and 4.2. We present a full proof in Section 4.3. The proof for the general case leads to a more effective 
implementation of Algorithm 1 than the proofs in Sections 4.1 and 4.2. 

Corollary 1. Given a finite presentation P for a 7L\E}\-module H, the word problem in *ZH is solvable. 

Proof. Given p = riihi + • • • + Tikhk- One can determine whether hi = hj since the word problem for H 
is solvable, and by the argument above, one can collapse p until it is no longer collapsible. □ 

We would like to remark that the 0(k 2 ) algorithm to collapse p = n\h\ + ■ ■ ■ + n^h^ is far from 
optimal. In putting the metabelian Fox matrix in upper-triangular form, many row operations may be 
required. This may well involve inverting many elements in Ki, and each time an inversion is done, we 
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must check that we are not dividing by zero. The expressions p = n\hi + ■ ■ • + n^h}~ get longer as more 
row operations are performed. An 0{k 2 ) operation done to many elements p of increasing length is very 
time consuming. Thus, a normal form for elements of H is preferred. A normal form for elements in a 
group is a canonical expression which allows for a quick decision of whether h = 1. For example, consider 
the polynomial multiplication (1 + 1 + 1 2 ) * (1 — t — t 2 ), and imagine for a moment that the terms it 1 ) are 
unrelated. Multiplying the polynomials using distributivity yields 1 * H — 1 * t — l*t 2 + ■ ■ ■ . One could 
then collapse the terms using the 0(k 2 ) algorithm described above. Now, remember that t 1 xt 1 = t 2 (i.e. 
t l is a normal form for a monomial). Then one may quickly multiply (l+t+t 2 )*(l—t—t 2 ) by grouping the 
terms in the product by degree, i.e. using the normal form: from the distributed product 1 — t — t 2 + ■ ■ ■ , 
one may rewrite the product (using a bit of memory) by scanning through the product only once. Here is 
the full product l-t-t 2 +t-t 2 -t 3 +t 2 -t 3 -t 4 = (1)1 + (-1 + l)t + (-1 - 1 + l)t 2 + (-1 - l)t 3 + (-l)t 4 . 
Assuming that integer addition is instantaneous, this method for collapsing polynomials is O(fc), where 
k is the length of the polynomial expression. 

4.1. The rational Alexander module. The classical integral Alexander module Aq(K) imbeds into 
the classical rational Alexander module Aq(K) by the map 

A%(K) = gW/gW G< 1 VG (2) ®iQ = Ao(K) 

which takes a i-> a <g> 1. This map is an imbedding since GW/G (2) is Z-torsion free [Cro63, Theo- 
rem 1.3]. The rational Alexander module Aq(K) is a module over Q a PID, and so A (K) = 
®"=i "2 i^ 1 ] I (Pi(t)) where pi{t)\pi + \{t) . Let d — Yl7=i deg(pi(t)). One can easily construct a Q \t ±1 ~\- 
module isomorphism 

n 
i=l 

by recording the polynomial coefficients in the factors. (There is an arithmetic formula for t-action on 
Q d which can be understood by the reducing of polynomials modulo the Pi(t) on the left-hand side). 

To sum up, there is an embedding of Z [t* 1 ] -modules G^/G^ ^ Q d , which induces an imbedding 
of rings 

Z [g (1) /G (2) ] ^Z[Q d ] 

The set Z[Q d ] is easy to work with. Its elements are functions from Q d to the integers (with finite 
support), which can be implemented on a computer using dictionaries keyed by <i-dimensional rational 
vectors and with values in the integers. 

Since recognizing € Q d is trivial, this method of tensoring H = G^/G^ 2 ' with Q gives a proof of 
Theorem 4 in the case G is the fundamental group of a knot compelement. 

From a computer programming point of view, this approach is untenable, insofar as computing 5\(K) 
is concerned. Encoding Z iG^ 1 ' /G^- 2 '] as a submodulc of Z[Q d ] is quite simple, but then arithmetic 
operations in Q d are required to upper-triangularize the metabelian Fox matrix. To avoid numerical 
round-off errors, one must use arbitrary precision arithmetic in Q d , which has substantial overhead. 
Arbitrary precision arithmetic uses rational numbers with arbitrarily large numerators and denominators. 
When a modern-day computer does an arithmetic operation an integer with greater than 20 digits 
(20 « log(2 64 )), it splits the integer into shorter pieces, which uses memory and takes time. 

We attempted implementing this 'tensoring with Q' approach in step 4 of Algorithm 1, and our 
program exceeded the computer's 8GB of RAM during computations for Si of several nine crossing 
knots. Perhaps one explanation is that even relatively few additions in Q can complicate the numerators 
and denominators: ^ — ^ + — ggggQ ■ The t-action also complicates the rational numbers (cf. 
Example 5). 

4.2. The Silver- Williams structure theorem. Given two maps / : U — > A, g : U — » B in the 

category of abelian groups, one can form the amalgamated sum 

A®uB:=A® B/{(f(u),0) = (0, g(u)) for all ueU] 



12 



PETER D. HORN 



The amalgamated sum is a pushout in the category of abelian groups. If both / and g are injective, 
then A and B imbed into A (Bu B. If one considers the case A = B, then one can iterate the operation. 
If A — B and / and g are injective, then this iteration is associative, i.e. 

{B ®uB)®uB-B ® v (B ®u B) 

One can form infinite amalgamated sums as well: • • • ®u B ®u B @u " ' ■ 

Given a finitely presented Z[Z]-module H which is Z-torsion free, Silver and Williams [SW12] construct 
explicit finitely generated abelian groups U and B and explicit maps f,g:U—¥B, By replacing U with 
U/Qserf + kerg) and B with B / (g(ker f) + /(kerg)) (and / and g by the induced maps) repeatedly 
until both / and g are injective. This procedure terminates after finitely many steps by the Noetherian 
property of finitely generated abelian groups. Silver and Williams construct a Z[Z]-module isomorphism 
®u B (Bu B (Bu ■ • • • The i-action on the amalgamated sum is given by shifting the B factor to 
the right. 

Using this isomorphism, we briefly describe an algorithm to decide whether a given clement of H is 
zero. Given x € H , apply the Silver- Williams isomorphism. The image lies in a finite range of summands 
C := B (Bu • • • ®u B. Since the amalgamating maps are injective, C <—> ■ ■ ■ (Bu B (Bu B (Bu ■ • ■ • Deciding 
whether an element of B (Bu 1 ■ ■ ®u B is zero is relatively simple. Let A — B (Bu ■ 1 ■ ®u B, which has one 
fewer B factor than C, so that C = A(Bu B. Then an element y £ C is zero if and only if it can be written 
as a sum of elements of the form (f(u), —g(u)), for u £ U. In other words y g C is zero if and only if the 
B-component can be pulled back by g and pushed into A via / so that the sum of the A-component with 
this 'swept back' i3-component is zero in A. One applies this sweeping (left) procedure until A = B, a 
finitely generated abelian group, where the detection of is simple. We summarize this in the following 
proposition, which with the above remarks gives a solution to the word problem for the Silver- Williams 
group presentation of H. 

Lemma 2. An element y of B (Bu 1 1 ' ®u B is zero if and only if y can be 'swept' to an element y' in 
the left-most B -factor and this y' = 0. 

While this approach allows one to detect in the ring Z [G«/G (2) ], it does not give a normal form 
for elements of G^ 1 ' /G^ 2 ' . For example, an element in C could be swept to the left out of C - there is no 
canonical B- factor to stop the sweeping. Thus, it does not speed up the 0(k 2 ) algorithm for collapsing 
elements of Z [G^/G^] , and so we chose not to implement the Silver- Williams method to complete 
step 4 of Algorithm 1 . 

4.3. Grobner bases. Grobner bases over Laurent polynomial rings have been used for algorithmic 
computations of the elementary Alexander ideals [GVHHUE06] but not, to the author's knowledge, to 
compute the entire Alexander module. 

Let R denote the polynomial ring Z[r, s]. We will use the degree reverse lexicographic order on 
monomials in R: r a s b < r c s d if and only a + b < c + d or {a + b = c + d and d < b }}. For example, 
s < r < s 3 < rs 2 < r 2 s < r 3 . This is a total ordering on the monomials in R so that any polynomial 
may be written in the unique way so that the monomials are increasing. A Grobner basis B for an ideal 
/ in R is a finite, list of polynomials fx , . . . /„ so that division of any polynomial R by the elements of B 
gives a unique remainder. Here, division of / by B means, roughly, to reduce / modulo the elements of 
B until it can be reduced no further. More specifically, division of / by B is achieved by running through 
this loop: while there is an /j g B such that / = q ■ /j + g for some q ^ and deg(g) < deg(/), set 
f = g and repeat until q = for all choices of i. The computations of the quotients q of course depends 
on the choice of fi for each pass through the loop, but the output of the loop, i.e. the remainder of / 
divided by B, is unique since B is a Grobner basis. We denote f%I the reduction of / € R modulo the 
Grobner basis B. It follows that f%I = if and only if / G I. 

One can construct an ordering on R d from a monomial ordering in R and a 'vector ordering.' Let 
v and be vectors in R d , and let i (respectively, j) denote the index of the last nonzero entry of v 
(respectively, We say v < ~w if either i < j or { i = j and v (i) < v$(i) }. This gives a total 
ordering on the elements of R d . Recall that the division algorithm in R computes quotients by looking 
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at the monomials in descending order. This division algorithm can be extended to R d by defining a 
monomial in R d to be an i?-monomial times a standard basis vector in R d . 

The notion of Grobner bases have been extended to finitely presented i?-modules. We refer the reader 
to [AL94, Chapter 3] but discuss the salient points here. Given a finite presentation matrix P for an 
i?-module M, one obtains an isomorphism M = coker P. Thus, the elements of M may be thought of 
as vectors in R d modulo the column space of P. A Grobner basis for the module M (with respect to the 
ordering on R d ) is a list B of vectors v\,...,v n such that the reduction of any v modulo B produces a 
unique remainder. We denote the reduction of v modulo the Grobner basis B by v%B. Such a basis 
is useful because it gives a normal form for elements of M, i.e. it = as elements of M if and only if 
lt%B = ~v$%B in R d . In particular, a Grobner basis affords one an easy method to detect € M, and 
such a basis is algorithmically computable. 

We are interested not in the ring R but in its quotient Z[r, s]/(rs — 1) = Z [t ±x ] ■ Given a presentation 
matrix P for a Z [t ] -module M, change all occurrences of t to r and t^ 1 to s. Then augment P by 
the matrix (rs — 1)1, and call the result P' . Reducing a vector it modulo a Grobner basis B' for the 
module presented by P' effectively reduces v by the columns of P and sets s = r^ 1 . In other words, 
for any v g R d , v%B' is a normal form for the image of v in M. We summarize this discussion in 
the following lemma. 

Lemma 3. Let H be a finitely presented Z [i ±:L ] -module. Then there exists a Grobner basis B for H 

that provides a normal form x%B € Z f^ 1 ]" 1 for an element x E H. Furthermore, x%B = if and 
only if x = G H . 

While polynomial division takes CPU time, this Grobner basis method has provided a more effective 
implementation of Algorithm l's step 4 than the 'tensoring with Q' method. While this method is 
slightly slower for the trefoil and figure-eight, it uses far less RAM and is able to calculate <5i for the 
knots 11„67, 12„293, and the (at most) 23 crossing knot from Example 4. We will illustrate the overhead 
benefits of the Grobner basis method. 

Example 5. Let P be the matrix 

^3 3i 2-t 



l-2t 3 - it j 

as in Example 4. Let us view P : Z[Z] 2 -> Z[Z] 2 by ~xt i Y ~atP. So the module P presents is Z[Z] 2 
modulo the row space of P. This agrees with the convention for Aq(K) as established in Example 3. 
We aim to find the normal form of the vectors x = (3t 3 + 2t 2 + 2)ei + (5t 2 — t i )e2 and t * x using the 
rational Alexander module and Grobner basis methods. 

Working over Q[Z], we calculate the Smith normal form of P to be D = LPR, where 

1 N \ „„a d- f 1 -7/3t + 5/3 



D= {o le-m + i) andi?= vo i 

For this example, L is irrelevant but may be computed using inverses. As a Q[Z]-module, M €5 Q is 
Q[Z]/(7t 2 - 13t + 7), though recall that M merely injects into M <E> Q as a Z[Z]-module. This quotient 
ring is isomorphic to Q 2 by recording the coefficients of the linear and constant terms of a polynomial's 
reduction modulo 7i 2 — 13t + 7. 

Now given any element y G Z[Z] 2 which projects to an element y € M, we have that <fi(V) — V^i where 
4> is the map that takes M to Q[Z]/(7t 2 — 13i+7) and ~ denotes taking the quotient modulo P or D. After 
composing with the isomorphism to Q 2 , we see that our it M> ( "147 ° ; llf) an< ^ ^x ^ ( "i V29 ' 
Notice how the i-action complicates the rational entries. Addition in this module complicates these 
rational numbers even more. 

Let r = t and s = t~ x . Using Macaulay2 [GS], we compute a Grobner basis as described above to be 

G -{{ )\ )>{ H 3 )\s+l)\ r + 1 
Macaulay2 computes lt%G = (-r 3 + 3r 2 -2r-14s + 14,0) and (rlt)%G = (-r 4 + 3r 3 -2r 2 -14s + 12, 0). 
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In the <g)Q method, the addition and t-action have inherent definitions in Q 2 . This is easy to implement 
and feels like it should be fast. However, the rational numbers become complicated after a few operations. 
Each subsequent operation takes increasingly more time, as evidenced by Example 6. In practice, this 
method's implementation exhausts our machine's memory during even simple <5i calculations. 

In the Grdbner basis method, the addition and i-action do not have inherent definitions, i.e. (x + 
f)%G ^ lt%G + f%G and (ff)%G ^ t(lt%G). After every addition and i-action, one must divide 
the result by the Grobner basis, which takes CPU time. However, the form ~x*%G is simpler - requiring 
less RAM than its Q 2 counterpart - and we find the Grobner basis implementation of encoding the ring 
Z [G«/G (2) ] more effective. 

Example 6. The i-action on Q 2 in the above example is given by t * (a, b) — (b + ^a, —a) . Consider 
the elements v = (-1420/147, 281/21) and x = (1282/147,-167/21) of Q 2 . These are the images under 
the isomorphism Z[Z] 2 /row P -> Q 2 of the elements (3t 3 + 2t 2 + 2,5t 2 -t 3 ) and (t 2 - 3i + 7, At 2 + 5t 3 - 2), 
respectively. Define a sequence v % by v n = t * v n -i + x. 

In the algorithm to compute Si, one must perform row operations on a matrix whose entries lie in 
Ki [t ±x ] ■ The coefficient of a monomial of one entry of this matrix is a formal quotient of elements of 
Z [GW/G' 2 )] . The numerator and denominator of this quotient may be quite long. Each term in the 
numerator is an integer times an element of G^/G^. If we think of G^/G^ as a subset of Q 2 , it is 
tempting to absorb this integer into the vector, but this is not permissible in the group ring Z[Q 2 ]. A 
single row operation on this matrix necessitates adding many monomials and hence many coefficients 
in Ki. Each addition a/b + c/d = (ad + bc)/bd in Ki involves three multiplications in the group ring 
Z[Q 2 ]. Each group ring multiplication involves many group operations in Q 2 , depending on the number 
of terms in the Z[Q 2 ]-elements which are multiplied. Note that each row operation on the matrix vastly 
increases the number of group operations in Q 2 that must be done. After a few row operations, there are 
a multitude of G^/G^ elements involved. If we think of these as elements of Q 2 , the numerators and 
denominators of the entries of each vector in Q 2 grow longer as the number of Q 2 operations increases. 
This, we believe, is why the 'tensoring with Q' implementation of step 4 exhausts our computer of all 
available memory. 

The sequence v n is intended to indicate how complicated the rational numbers become after many 
arithmetic operations. We do not foresee the v n arising in a Si computation. While the rational model 
for G^/G^ offers a faster implementation for the Z [G^/G^j-operations than the Grobner basis 
model - at least when few operations are required - the rational model slows down significantly when 
many operations are performed. We computed v n to n — 50000 and noticed a significant slowing down 
of the computation as n increased. See the table in Figure 6. The CPU time required to compute v n 
appears to be exponential in n. 

Note that the terms of this sequence v n written in the Grobner basis model will also become com- 
plicated as n increases. By the Grobner basis in Example 5, the second entry in v n will be an integer, 
and the first entry will be a polynomial in r (with perhaps one or two s terms) whose degree grows sub- 
linearly in n and most (all but two) of whose coefficents have absolute value less than 7/2 (Macaulay2 
computes a quotient so that the remainder's leading coefficient is as small in absolute value as possible, 
i.e. 13r % 5r — — 2r). Note that the second entry has absolute value less than 3/2. Thus, v n can be 
stored in a computer by approximately n digits. This is more efficient than the rational method, at least 
experimentally, for approximately n digits were used just for the denominator of one entry in the Q 2 
representation of v n , according to the table above. 

5. Computations 

Using a 2.5 GhZ Intel i5 processor with four cores and 8 GB of RAM, we made the following Si 
calculations. The implementation is discussed in Section 6. We obtained the PD code for all knots using 
Knotlnfo [CL13]. For the 23 crossing Example 14, we used the program Link Sketcher, which we found 
on Knotlnfo, to compute the PD Code. 
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0j . CPU time to compute 

btep i . , 

^loooi; ln seconds 

1 0.048321 

2 0.146888 

3 0.315882 

4 0.57017 

5 0.921086 
10 4.870255 
20 28.885591 
30 84.075005 
40 181.961961 

45 249.442993 

46 265.189638 

47 281.562427 

48 298.729361 

49 316.309955 

50 334.467209 



Digits in denominator of 
first entry of «ioqqi 

844 

1689 

2534 

3380 

4225 

8450 
16901 
25352 
33803 
38029 
38874 
39719 
40564 
41409 
42254 



Example 7. The trefoil knot 3i has PD code X[l, 5, 2, 4], X[3, 1, 4, 6], A[5, 3, 6, 2]]. From this PD code, 
our program calculated <5i(3i) = 1 in 0.2128 seconds. 

Example 8. The figure-eight knot 4i has PD code X[4, 2, 5, 1],X[8, 6, 1, 5],X[<o, 3, 7, 4], X[2, 7, 3, 8]. Our 
program calculated <5i(4i) = 1 in 0.5616 seconds. 

Example 9. The 2, 5-torus knot 5i has PD code X[2, 8, 3, 7], X[4, 10, 5, 9], X[6, 2, 7, 1],X[8, 4, 9, 3], [10, 6, 1, 5]. 
Our program calculated #i(5i) = 3 in 9.5816 seconds. 

Example 10. The knot 5 2 has PD code X[l, 5, 2, 4], X[3, 9, 4, 8], X[5, 1, 6, 10], X[7, 3, 8, 2], X[9, 7, 10, 6]. 
Our program calculated <5i(52) = 1 in 1.4597 seconds. 

Example 11. The stevedore knot 6 1 has PD code X[l, 7, 2, 6], X[3, 10, 4, 11], X[5, 3, 6, 2], X[7, 1, 8, 12], 
X[9,4,10,5],X[11,9, 12,8]. Our program calculated <5i(6i) = 1 in 2.7247 seconds. 

Example 12. The knot 11„ 67 has PD code X[4,2,5, 1],X[8,4,9,3], X[ll, 17, 12, 16],X[14,5, 15,6], 
X[6, 15, 7, 16], X[9, 19, 10, 18], A[17, 11, 18, 10], [19, 1, 20, 22], X[13, 20, 14, 21], X[21, 12, 22, 13], X[2, 8, 3, 7]. 
This PD code produces a Wirtingcr presentation with 11 generators, and our program would attempt 
to put a 10 x 10 matrix over Ki [t^ 1 ] in upper-triangular form. However, ll n 67 is a 3-bridge knot, 
and thus its fundamental group has a presentation with 3 generators and 2 relations. Using Tietze 
transformations in GAP [GAP12], we were able to reduce the Wirtinger presentation to the gen- 
erators x±,X2,X3 and relations x^ x%x^ x^ XiX^x^x^ Xtf^ x~^ XiX~^ x§~ x^[ x\x^x^ XiX^x^xl^ and 
X\X^ x^ x^x^ xl^ x^ x^ x\x%x^ x\x%x\ x^x^ x^ Xix^x^x^x^ . In this presentation, each Xi maps 
to a generator of Z, and so a presentation with a 'nice' (as in Example 3) generating set can easily be 
obtained. This results in a 2 x 2 metabelian Fox matrix that must be put in upper triangular form. Of 
course, the entries in the 2x2 matrix are more complicated than in the original 10 x 10 matrix, but the 
5i computation is faster. 

Given the reduced Wirtinger presentation, our program calculated <5i(ll n 67) = 3 in 11324 seconds. 
As mentioned around Theorem 2, this is the first known knot in the tables with Si > Sq. 

Example 13. The knot 12„ 293 has PD code X[l, 4, 2, 5], X[3, 10, 4, 11], A[5, 12, 6, 13], A[16, 8, 17, 7], 
X[9, 2, 10, 3], X[ll, 8, 12, 9], X[20, 13, 21, 14], X[(a, 16, 7, 15], X[24, 17, 1, 18],X[22, 19, 23, 20], X[U, 21, 15, 22], 
X [18, 23, 19, 24]. Using GAP, we found a simplified presentation for the fundamental group of the 
knot complement with generators a, b, c and relators b 2 ab~ 2 a~ 1 c~ 1 bacab 2 a~ 1 b~ 2 a~ 1 c~ 1 a~ 1 b~ 1 cab and 
b~ a -1 c~ 1 bab~ 1 cabcab^ 1 a -1 c~ 1 bab~ 1 a -1 c~ 1 ba~ 1 b~ 1 caba^ 1 cr x . Here a maps to a generator under abelian- 
ization, and b and c map to zero. Our program calculated <5i(12„293) = 3 in 8393 seconds. This is another 
example of a knot with 2 • genus = 1 + Si > 2 ■ Sq . 
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Example 14. Let K denote the knot from Example 4. The diagram for K in Figure 1 has 23 crossings, 
and its PD code is X[M, 1, 35, 2], X[2, 33, 3, 34], Jf [32, 3, 33, 4], X[4, 31, 5, 32], X[30, 5, 31, 6], X[6, 29, 7, 30], 
X[7, 40, 8, 41], X[8, 20, 9, 19], X[22, 10, 23, 9], X[37, 10, 38, 11],X[11, 36, 12, 3,7], X [12, 24, 13, 23], 
X[46, 13, 1, 14], X[U, 45, 15, 46],X[44, 15, 45, 16], X[16, 43, 17, 44], X[42, 17, 43, 18], X[18, 41, 19, 42], 
X[27, 20, 28, 21], X[21, 26, 22, 27], A[24, 36, 25, 35], A[38, 26, 39, 25], X[28, 40, 29, 39]. Using Tietze trans- 
formations in GAP, we were able to reduce the Wirtinger presentation with three generators X\,x^,x^ 
and three relations x\x^ 1 x\X2X^ 1 X2XJ 1 X2X^ x^ 1 x^x^ 1 and 
Our program calculated 8\(K) = 1 in 0.7698 seconds. 

6. An implementation of the algorithm 

We implemented Algorithm 1 in Sage 5.3. Most of the classes and functions are written in Python 
2.7.2, and we used Macaulay2 version 1.4 to compute the Grobner basis for Aq(K). Sage incorpo- 
rates Python and Macaulay2, and so our program can be executed in Sage with one click of a button. 
Sage [S + 12], Python [Pytl2], and Macaulay2 [GS] are open source and cross-platform. 

Our implementation can take as an input the PD code for a knot diagram or a Wirtinger presentation 
for the knot complement's fundamental group. We use a script written by REU students to convert 
the PD code to a Wirtinger presentation [REU11]. At this point, the algorithm runs quite similarly to 
Example 3. After switching to a convenient generating set for the fundamental group, the Fox matrix is 
then calculated, and a presentation matrix (the abelian Fox Matrix ) for GW/G (2) £ Al(K) is recorded. 
The abelian Fox matrix is fed to Macaulay2, which computes a Grobner basis for the classical Alexander 
module. Separately, the entries of the original Fox matrix are written modulo G^ 2 \ using a splitting 
which treats the first generator of the Wirtinger presentation as t, so G/G^ = G^/G^ x Z, as in 
equation 2. Our program finds a column to delete by reducing the G^/G^ -generators modulo the 
Grobner basis computed earlier, and deletes it, yielding the metabelian Fox matrix, which presents 
A\(K) as a Ki [i* 1 ] -module. Finally, this metabelian Fox matrix is put in upper-triangular form by 
row and column operations, and the degrees of the polynomials are computed and added, giving 5\{K). 

We wrote classes for the elements of -k\, G/G^ 2 \ ZTi, Ki, and Ki [t^ 1 ] , and functions to add and/or 
multiply elements in each of those sets, which culminated in the polynomial ring operations. Throughout 
the upper-triangularization process, we need the degree function in Ki [t^ 1 ] , which amounts to recog- 
nizing when the coefficient of a term is zero or not, hence the importance of Section 4. In Ki [t* 1 ], 
we need to compute quotients of polynomials on the left and right so that we can do row and column 
operations to upper-triangularize the metabelian Fox matrix. 

We implemented group ring elements in Z [G^/G^] as Python dictionaries. A dictionary is a 
function whose domain consists of keys. The range of the function is the set of values for the dictionary. 
In Python, keys must be hashable, so that dictionary lookup time (recalling the value of the function 
given a key) is constant time. A group ring element in Z [GW/G (2) ] can then be thought of as a 
dictionary keyed by elements of G^/G^ with values in Z. Multiplication of elements in Z [GW /G^] 
follows the rule (ax + by)(cz + dw) = ab(x + z) + ad(x + w) + bc(y + z) + bd(y + w) (we use here that 
the group G^'/G*- 2 -* is abelian). Python dictionaries can be used to multiply and group like terms, 
since we have a normal form for elements of G^/G*- 2 ). We include the pseudo-code for our group ring 
multiplication function: 

def multiply_dictionaries (dictl , dict2) : 
product = O ; 
for v in keys of dictl: 

for w in keys of dict2: 

# v and w are vectors reduced by the Groebner basis GB; 
product_key = (v+w) % GB; 
product [product_key] = 0; 
for v in keys of dictl: 

for w in keys of dict2: 



ON COMPUTING HIGHER-ORDER ALEXANDER MODULES OF KNOTS 



17 



product [(v+w) '/. GB] = product [(v+w) '/, GB] + dictl [v] *dict2 [w] ; 
for k in keys of product: 
if product [k] == 0: 

delete product [k] ; 
return product 

The first for loop creates all of the sums of group elements that occur in the product in the group 
ring. The second for loop computes all the integer coefficients in the product. The third for loop removes 
all terms in the product with coefficient zero. If I is the length of the first factor and r is the length of 
the second factor, this function runs through Zlr loops, at most. The returned product is completely 
collapsed. Alternatively, one could implement group ring elements as lists. Multiplication and collapsing 
of lists must run through at most Ir + (lr) 2 loops. 

A Sage worksheet for computing S\ is available at https : //pdhorn. expressions . syr . edu/sof tware/. 
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